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Abstract 



The sources of W + n-jet events in hadron collisions are higher-order QCD 
processes, but also multiple-parton interactions. A subprocess producing a 
W + /c-jet final state, followed by one producing I jets in the same nucleon- 
nucleon interaction, will result in a W+n-jet event if k+l = n. In the simplest 
case a PF + 2-jet event can be produced by a quark-antiquark annihilation into 
W and a 2-jet event occurring in the same proton- antiproton interaction. We 
compute that this happens at the 10% level of the higher-order QCD processes 
for the type of cuts made by the Tevatron experiments. For jet pt values of 
order 5~10 GeV, multiple-parton interactions dominate higher-order QCD- 
processes. The emergence of this new source of W + n-jet events towards lower 
Pt simulates the running of a s ; it is imperative to remove these processes 
from the event sample in order to extract information on the strong coupling 
constant. Also, BFKL studies of low-p? jet cross sections are held hostage to 
a detailed understanding of the multiple-parton interactions. We perform the 
calculations required to achieve these goals. A detailed experimental analysis 
of the data may, for the first time, determine the effective areas occupied by 
quarks and gluons in the nucleon. These are not necessarily identical. We 
also compute the multiple-parton contribution to Z + n-jet events. 



I. INTRODUCTION 



The very clean samples of W, Z + n-jet events produced in proton-antiproton interactions 
at the Tevatron |IJ have become a laboratory to study QCD ||. In principle, these events 
can be used to determine the strong coupling a s . In this paper we point out that the same 
events can be used to study multiple-parton interactions. It would, in fact, be treacherous 
to determine a s without removing their contribution to the data sample. We will show 
that their contribution strongly increases relative to the higher-order QCD processes with 
decreasing values of pt- This behavior obviously mimics the running of the strong coupling. 
Furthermore, studies of BFKL physics with jets cover the low-pr region where multiple- 
parton processes eventually dominate and will have to await a detailed understanding of 
these processes. 

We have computed the contributions of multiple interactions and higher-order QCD 
processes to the processes W, Z + 2-jets; see Fig. [I]. In the simplest case of a W + 2-jet 
event the final state can be produced by a quark-antiquark annihilation into W and a 2-jet 
event, both occurring in a single beam crossing ||. The two-parton subprocess may involve 
the same, or different beam particles. In most cases the latter process can be identified by 
the detector and removed from the data set H. Our conclusions are that: i) higher-order 
QCD processes and multiple-parton interactions are similar in magnitude for values of 
the jets of order 5~10 GeV and ii) even for jets with a minimum px of 20 GeV the multi- 
parton contribution is at the 10~20% level. Notice that a contribution of this magnitude is 
sufficiently important to interfere with a determination of a s . The obvious way to remove 
it is to study the event sample as a function of the minimum p? of the jets. We have also 
calculated the multiple interaction and higher-order QCD contributions to W + 3-jet event 
and comment on their relative importance when the number of jets produced in association 
with the W increases. Finally, we repeated the calculations for Z + 2, 3-jet events. 

Multi-parton scattering is of theoretical interest because it probes the partonic structure 
of hadrons in novel ways. The joint probability that a quark and antiquark annihilate into a 
W, and that another pair of partons produces a 2-jet final state in a single proton-antiproton 
interaction is 

aw 2-jet = crw • i 1 ) 

The result is obvious except for the factor — It takes into account the fact that, once the 
first parton process has taken place, the area of the beam has been reduced to the area of the 
proton in which the first interaction occurred. In principle, there could exist correlations 
in the double-parton-scattering structure functions that would invalidate this expression. 
For instance, we are neglecting the correlation due to longitudinal momentum conservation, 
which should introduce factors of the form (1 — x\ — X2), with Xi being the fraction of 
momentum carried by the parton i. In the small x regions, covered in our calculations, 
Eq. (HD is valid. The effective area ttR 2 should be of order of the inelastic cross section, 
or about 50 mb. Earlier studies of final states with multiple jets and multiple Drell-Yan 
pairs [0 , although admittedly plagued by difficult systematics or low statistics, have agreed 
on the much smaller value of roughly 15 mb. This may indicate that hard quarks and gluons 
only occupy a small area of the high-energy proton. Anticipating the higher quality of the 
Tevatron data, it may ultimately be possibly to separately measure the areas occupied by 
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quarks and gluons. Several authors have speculated that gluons may be clustered in the 
vicinity of valence quarks 0. 



II. W + 2-JETS EVENT RATES 

The multiple-interaction cross section for the production of a final state A in association 
with the state B, is given by 



a A -B = °a°b 



1 , 2 / 1 (N-l) | (N-l) 



TtR 2 V KR 2 (Tinel <7 inc l 



(2) 



for A / 5. Here iV stands for the average number of interactions per beam crossing. 
Summation over A and B states leading to the same A — B final state is understood. The 
first term corresponds to double-parton interactions involving the same beam particles, like 
in Eq. ([J), while the last term describes simultaneous interactions between different beam 
particles. As for the second term in Eq. (1), it corresponds to the interaction of two partons 
in a proton in one bunch with two different protons in the other bunch. It can be understood 
as the interference term for the first two reaction mechanisms. In our calculation we use 
R = 0.7 fm, which corresponds to an effective cross section ttR 2 = 15.4 mb, cr ine i = 51.7 mb, 
and N = 1.8. This value of R is also consistent with the the one derived from the analysis 
of the 7 + 3-jet production by CDF [f|. 

In order to obtain the multiple-parton and QCD cross sections for the production of a W 
and two jets, we used the package MADGRAPH [[F] to evaluate the relevant tree-level matrix 
elements. These agree with previous evaluations. We imposed a minimum-pr cut on the 
transverse momenta of the jets and required that the jets (partons) be in the rapidity region 
\r]j\ < 3, and separated by AR = ^At] 2 + A~cj? > 0.9. We used the CTEQ3M structure 
functions || and evaluated the QCD scales at s/2, where § is the subprocess center-of-mass 
energy. 

Our results for the W + 2 jets production are exhibited in Fig. [2] as a function of the 
minimum-pr cut. We present in Fig. the pr spectrum of the jets produced in association 
with the W for both multiple-interaction and higher-order QCD processes. As can be seen 
from this figure, the jets produced in multi-parton scattering are much softer than those 
produced by higher-order processes. This is why multiple-particle interactions become more 
important with decreasing p?. 

In Fig. f| we plot the statistical significance (signal over square root of background) for 
multi-parton interactions assuming an integrated luminosity of 100 pb _1 , Br(W / — ► ev e ) = 
10.8%, and an efficiency of 34% to detect the W [l|]. The double-parton interactions can 
be observed in the W + 2-jet final state with a significance in excess of 5a provided the 
minimum pt of the jets is smaller than 10-12 GeV. Clearly, the signal for double-parton 
scattering can be further enhanced, e.g. by a cut requiring that the jets are back-to-back in 
the transverse plane in a W + 2 jet event. The poorly known effective area ^2 dominates 
the uncertainty of the calculations. The results shown assumed a common value of 15.4 mb 
for quarks and gluons, corresponding to R = 0.7 Fermi. It is clear from these results that 
the Tevatron has the opportunity to determine this quantity with much improved precision. 
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III. W + 3-JETS EVENT RATES 



Double-parton interactions can also be studied in the production of W + 3-jets. The 
cross section for this process is given by Eq. (0) with 

The cross sections for this final state via higher-order QCD, double-parton interactions and 
multi-beam interactions are shown in Fig. [5] as a function of the minimum p? of the jets. As 
we can see from Fig. [], the transverse momenta of the W produced via higher-order QCD 
are larger than those generated in the double-parton interaction. This suggests that a cut 
like p?f < 20 GeV should increase the signal-to-background ratio. 

Figure (|also shows the statistical significance for the double-parton production ofW + 3- 
jets as a function of the minimum transverse momenta of the jets, indicating that the signal 
should be visible for p™* < 12 GeV. This curve was obtained using the same luminos- 
ity and efficiencies assumed for the W + 2-jet signal. No further cuts were introduced to 
enhance the signal. Notice that double-parton scattering is observable in both W + 2-jet 
and W + 3-jet channels for approximately the same range of transverse momenta. There- 
fore, the observation of a signal in one channel can be confirmed in the other. Moreover, 
it is possible to observe the presence of double-parton scattering by analyzing the ratio 
R3/2 = 0"VK+3-jets/cTy+2-jets- As shown in Fig. double-parton interactions enhance this ratio 
by a factor 1.5-2 for \ow-pr cuts. 



IV. Z + n-JETS 

Multiple-parton scattering can also be studied through the production of a Z accompa- 
nied by 2 or 3 jets. In order to study these processes we must evaluate 

^ACr B \z2-ict =VzV2-jet , (3) 
^filzs-jet = CTZjetO-2-jet + CTzO-3-jet • (4) 

We exhibit in Table [I] ([TT]) the cross sections for the production of Z + 2-jets {Z + 3-jets) 
from higher-order QCD processes, double-parton scattering and multiple-beam interactions. 
These results were obtained imposing a minimum-p-^ cut on the transverse momenta of the 
jets. As before, we required that the jets (partons) populate the rapidity region \r\j\ < 3, 
separated by AR > 0.9. 

The statistical significance of the multi-parton contribution to the Z + 2-jet production 
is shown in Fig. |8] for Z decaying into e + e _ or pairs, and an integrated luminosity 

of 100 pb _1 . The Z acceptance was taken to be 41% |fj]. We conclude that the multi- 
parton scattering can be observed at a 3a level for p^ ut < 14 GeV. The signal should be 
very clean for low values of the px cut, even without further cuts requiring the Z and jets 
to have pair-wise balancing transverse momenta. The signal for multi-parton scattering in 
the production of a Z + 3-jets is also visible for a px cut smaller than 13 GeV, as shown 
in Fig. |8|. Analogously to the W + ra-jets case, the signal for multi-parton interactions in 
Z + n-jets also appears in the ratio R 3 / 2 = crz+3-jcts/ cr z+2-jets, as shown in Fig. ||. 
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TABLES 

TABLE I. Z + 2-jet cross sections at the Fermilab Tevatron, with the cuts and parameters 
discussed in the text. 



p T cut (GeV) 




^multMnt. (pb) 


mult. part, r i \ 
CJ Z2,jet (P b ) 


5 


800. 


1578. 


710. 


10 


269. 


159. 


71. 


15 


126. 


36. 


16. 


20 


67. 


12. 


5.3 


25 


40. 


4.7 


2.1 



TABLE II. Z + 3-jet cross sections at the Fermilab Tevatron, with the cuts and parameters 
discussed in the text. 



Pt cut (GeV) 


»SS (pi-) 


O-Z 3-jet (Pb) 


^mult^int. (pb) 


mult. part. / i \ 
^3-ict (P b ) 


5 


2582. 


195. 


906. 


408. 


10 


1386. 


42. 


49. 


22. 


15 


868. 


14.5 


7.2 


3.2 


20 


586. 


6.2 


1.6 


0.72 


25 


412. 


3.0 


0.46 


0.21 
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FIGURES 




(a) (b) (c) 

FIG. 1. W,Z + 2-jet multiple-parton scattering. 
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FIG. 2. Total cross section for the production of a W accompanied by 2 jets, as a function 
of the minimum-pr cut, for: QCD processes (open circles), double-parton interactions (triangles), 
multiple interactions (squares). 
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FIG. 3. pt distribution of the jets in mb/GeV, for: (a) QCD 2-jets production, (b) VF + 2-jets 
production. 
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FIG. 4. Statistical significance (Signal/ ^/background) of the double-parton interactions for the 
production of a W + 2-jets (squares) and W + 3-jets (triangles) as a function of the minimum 
transverse momenta of the jets. 
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FIG. 5. Total cross section for the production of a W accompanied by 3 jets, as a function of 
the minimum-px cut, for: higher-order QCD processes (open circles), double-parton interactions 
(triangles), multiple interactions (squares). 
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FIG. 6. pt distribution of the W in mb/GeV, for: (a) W + 1-jet production; (b) QCD 
W + 3-jets production with the solid line standing for the total while the dashed, dotted, and 
dot-dashed lines stand for the contribution from the quark-gluon, quark-quark, and gluon-gluon 
subprocesses respectively . 
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FIG. 7. #3/2 in H^ + n-jet production as a function of the jet minimum transverse momentum. 
The squares stand for the sum of the QCD higher-order and double-parton contributions, while 
the triangles represent only the higher-order QCD contributions. 
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FIG. 8. Statistical significance (Signal/ ^/background) of the double-parton interactions for the 
production of a Z + 2-jets (squares) and Z + 3-jets (triangles) as a function of the minimum 
transverse momenta of the jets. 



13 



0.5 



0.45 - 



0.4 - 



0.35 - 



0.3 - 



0.25 - 



0.2 



0.15 - 



0.1 



0.05 - 




10 



15 



20 



25 30 
Pt (GeV) 



FIG. 9. -R3/2 i n Z + n-jet production as a function of the jet minimum transverse momentum. 
The squares stand for the sum of the QCD higher-order and double-parton contributions, while 
the triangles represent only the higher-order QCD contributions. 
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